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Abstract 

Ratios of multiplicity moments, Hq (cumulant over factorial mo- 
ments Kq/Fq), have been observed to show an oscillatory behaviour 
with respect to order, q. Recent studies of e^e~ annihilations at LEP 
have shown, moreover, that the amplitude and oscillation length vary 
strongly with the jet resolution parameter ycut- We study the pre- 
dictions of the perturbative QCD parton cascade assuming low non- 
perturbative cut-off {Qq ~ Aqcd ~ few 100 MeV) and derive the 
expectations as a function of the cms energy and jet resolution from 
threshold to very high energies. We consider numerical solutions of 
the evolution equations of gluodynamics in Double Logarithmic and 
Modified Leading Logarithmic Approximations (DLA, MLLA), as well 
as results from a parton MC with readjusted parameters. The main 
characteristics are obtained in MLLA, while a more numerically accu- 
rate description is obtained by the MC model. A unified description of 
correlations between hadrons and correlations between jets emerges, 
in particular for the transition region of small ycut' 



1 Introduction 



The multiplicity of hadrons is the simplest global characteristic of the final 
state in particle collisions. It was an early conjecture that the main trends of 
the mean multiplicity and higher multiplicity moments could be reproduced 
by the corresponding perturbative QCD calculations for the quark-gluon cas- 
cade. First results have been derived for the asymptotic high energy limit in 
application of the Double Logarithmic Approximation (DLA) |ip2p3j . Taking 
into account the leading contributions from the collinear and soft radiation 
singularities and angular ordering |4|5j one arrives at a probabilistic descrip- 
tion of the parton cascade which evolves from the high energy scale Q down 
to the hadronization scale Qq. The perturbative expansion in the coupling 
as can be resummed and the exponent can be expanded into a power series 
of ^/oTs. In this way results on mean multiplicity |5pHp(ij and higher moments 
|2pHp7j have been obtained, subsequently the next-to-leading log corrections 
in MLLA |8l9llUllI] and yet higher order terms |12pi3pi4] have been derived. 
All these results can be obtained as different approximations to the MLLA 
evolution equations for the generating function of quark and gluon jets [T5] . 
A full solution of these equations can be obtained by numerical methods 

cnini. 

One motivation of such studies was - and still is - to find out how far 
one can extend perturbative calculations into the regime of real strong in- 
teractions which determine multiparticle production. In this way, we hope 
to learn more about the colour confinement mechanism. The MLLA results 
on mean multiplicities and single particle spectra give a surprisingly good 
description of the data which led to the notion of "Local Parton Hadron Du- 
ality" (LPHD [9 ). The results here depend on only two parameters, the QCD 
scale A and a single non-perturbative parameter, the transverse momentum 
cut-off Qo with Qq > A; an overall normalization factor K is also allowed for. 
Such parton level calculations have been applied to a large variety of inclusive 
observables, in particular to multiparticle correlations ( "Generalized LPHD" 
|18j): even quasi-exclusive processes are within reach In general, one 

observes that with increasing accuracy, the perturbatively calculated observ- 
ables are in better agreement with the data. This simple model does not 
contain any hadronic resonances, so it is clear that it can only describe suffi- 
ciently inclusive observables and should not be considered for too restricted 
regions of phase space (for pertinent reviews, see |2UI21l22l23j ). 

Whereas there is no problem with inclusive single particle results, the 
status of higher multiplicity moments (factorial moments Fq and cumulant 
moments Kg, see below), which are integrals over the respective q particle 
correlation functions, is more controversial. The data on quark and gluon 



1 



jets for g < 5 P3] deviate with increasing order q from the higher order 
logarithmic approximations |25j; they correspond to an expansion in q^fal 
rather than ^Jal |2nf22j and the known terms fit the mean multiphcity (mo- 
ment g = 1) but the moments of higher order become much larger than 
the data. On the other hand, numerical solutions of the MLLA evolution 
equation yield very good agreement with these data |16|17j : remarkably, a 
common fit of hadron and jet multiplicities is possible and results in 

Ktot ~ 1 (1) 

for a small cut-off Qo > A where Ki^t is the ratio of the total hadron to parton 
multiplicities. This ratio Ktot ~ 1 was also found to be consistent with the 
quark and gluon multiplicities in high p-r jets observed at the TEVATRON 
|2f)j where the higher order MLLA corrections are taken into account as well. 

An interesting prediction from the asymptotic ^Jo^s expansion con- 
cerns the ratio of moments 

a, = KJF, (2) 
These moments show an oscillatory behaviour with the first minimum at 

g™n^T— + - + 0(7o), h, = - 7o = ^^ (3) 
/ii7o 2 24 TT 

where Nc = 3 and 70 denotes the leading order multiplicity anomalous di- 
mension at the considered energy scale, numerically g^m ~ 5(±1) at LEP 
energies. Such oscillations have indeed been observed in e~^e~ annihilations 
at SLC |2Zj and recently at LEP |28p29j with the first minimum at qmin ~ 5 
as expected. 

Since the absolute size of the moments are not in quantitative agreement 
with these high energy predictions the physical origin of these oscillations re- 
mains controversial. It has been noted that the truncation of the multiplicity 
distribution at large multiplicity, n, also gives rise to oscillations |3nj- Fur- 
thermore, a superposition of two components (like 2 and 3 jet events), each 
one modeled without oscillations, can lead to oscillations in the full sample 

m 

A new element has been brought into the discussion recently through 
analyses of the L3 measurements |28f29j of multiplicity moments for hadrons 
and narrow jets at high resolution (small t/cut parameter). The measure- 
ments feature a rapid variation of both the oscillation amplitude and oscil- 
lation length in q as function of i/cut- This behavior does not appear in the 
asymptotic solutions of the evolution equation [13]. 



2 



In this paper we aim at an understanding, within perturbative QCD, of 
the Hq oscillation phenomena in both hadron and jet final states. We study 
first the solutions of the evolution equations (both DLA and MLLA) with 
jet hardness from threshold up to the asymptotic regime; furthermore, we 
investigate the variations of the jet multiplicity moments with resolution, 
from low resolution (fat jets) up to high resolution (narrow jets) - ultimately 
up to the scale for the parton-hadron transition. These calculations depend 
only on the QCD scale. A, and on one parameter, Qc-, which refers to the 
transverse momentum cut-off. This transverse momentum cut-off is given by 
either the scale of the arbitrarily chosen jet resolution or the characteristic 
hadronic scale Qc = Qo- The full results including the dependence on Qc, 
are obtained by numerical methods. The need for an overall normalization 
factor K is considered for the mean multiplicity. 

A similar theoretical scheme is realized in a parton level MC (ARIADNE 
jH2!) which is based on sequential parton radiation from colour dipoles jH3p34j 
with a kr cut-off. We have readjusted the parameters A and Qq in this MC 
to describe hadronic final states without an additional hadronization phase, 
assuming again a duality between hadron and parton final states at scale Qq. 

2 Definition of moments and jet resolution 

The distribution, P„, of the multiplicity, ra, of particles or jets in an event 
can be characterized by its moments. One considers the factorial moments 
fq or the normalized moments Fq 



fq = Y,<n-l)...{n-q + l)P^, Fq = fq/N\ N = f, (4) 



with mean multiplicity A^. Furthermore, one introduces the cumulant mo- 
ments kg and Kg which are used to measure the genuine correlations without 
uncorrelated background in a multiparticle sample 



in particular K2 = F2 — 1, K3 = F3 — 3F2 + 2; for a Poisson distribution 
Ki = l, Kg = for g > 1. 

These moments can also be obtained from the generating function of the 
multiplicity distribution 
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by differentiation 

where Y denotes here some Icinematic variable hlce the total energy. 

In a high energy collision, many hadrons are produced and they are found 
to cluster typically into jets of particles. The number of such jets depends on 
the resolution, which is defined through an algorithm in terms of a resolution 
parameter. In this paper, we choose the so-called Durham algorithm [HK] . 
At resolution i/cut one combines particles into jets iteratively until all pairs 
of jets or particles have relative transverse momentum > ycutQ"^ for small 
relative angles 9jj, or, in general, 

y,, = 2min(E2,£;|)(l-cose,,)/Q2 > y,^f (8) 

In our application to e^e~ annihilation, Q is the total cms energy. At very 
low resolution {ucut 1) one sees only two jets which evolve from the primary 
qq system, whereas with increasing resolution more and more jets are resolved 
until at very high resolution {jjcut 0) all final state hadrons are resolved, 
so 

Ucut 1 : Njet = 2, Vcut : Nj^t = Nhadron- (9) 

In this paper we study both the multiplicity of jets and the multiplicity of 
hadrons and the respective moments. In particular, we are interested in the 
variation of these observables with energy and resolution as has recently been 
reported in the experimental works . 2 8, 2.9 j . 

In the popular Monte Carlo models of today, the parton cascade is ter- 
minated at a scale of Qo ~ 1 GeV; this perturbative phase is followed by 
the hadronization phase for which non-perturbative models, typically with 
a rather large number of parameters, are introduced. In our calculation we 
follow the idea of LPHD 9J and let the parton cascade evolve further down 
to a lower scale Qq of a few 100 MeV of the order of the QCD scale A itself. 
Then we compare the result directly to the hadron final state. In this way, 
besides the common parameter A only one extra non-perturbative parameter 
Qq is introduced which can be interpreted either as a typical hadron mass 
or as inverse hadron radius. In the theoretical calculations we obtain the 
hadronic final state for Q^ — > Qq and we replace Q by 

Ucut 1 : Njet = 2, Vcut {Qo/Qf ■ Njet = Nhadron- (10) 

An important feature is the transition from jets to hadrons for small 
resolution parameter tjcut which tests the hadronization models in the region 
of small momentum transfers, the so-called "soft" region. In this paper the 
consequences of such an approach for multiparticle correlations in the soft 
region is investigated in the theoretical models and compared with the data. 
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3 Solution of evolution equations for multi- 
plicity moments 

3.1 Double Logarithmic Approximation 

In this approximation we consider a single jet of particles or sub-jets at 
transverse momentum cut-off kj- > Qc] for e~^e~ annihilation, two jets have 
to be superimposed, with Qc then corresponding to the resolution i/cut- Only 
the most singular terms of the parton splitting functions are kept and recoil 
effects are neglected. Then we can restrict ourselves to gluodynamics, i.e. 
we neglect the production of quark anti-quark pairs which would represent 
a nonleading contribution. A common simplification of these calculation is 
the restriction to 1-loop results; furthermore, in this application we consider 
only light quarks. In this case, the generating function follows an evolution 
equation 



^(^-""Vo ^^'^o(y')[^(y',^)-l] (11) 
Z(0,m) = m (12) 

with the multiplicity anomalous dimension 

TT y + Xc 

y, = ln(i?/Q,), A, = ln(Q,/A), = b=^Nc-'^nf (14) 

b 3 3 

where N^. = 3. The argument of is related to the transverse momentum, 
i.e. y = ln{kT/Qc)- In our analytic calculations we take 

A = 300 MeV, A = ln(go/A) = 0.015, Uf = 3. (15) 

The initial condition ()12|) implies that at threshold {Yc = 0) the jet Z{Yc,u) 
contains only one particle (sub-jet). 

Differentiating the evolution equation (jllj) with respect to u at u = 1, 
one obtains the following equation for the multiplicity, N: 

N'{Y,)= r dy'^l{y')N{y') (16) 
Jo 

or, equivalently, 

N"{Y,) - ^^,{Y,)N{Y,) = 0, N{0) = 1, iV'(O) = 0. (17) 
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Figure 1: Multiplicity, A^, in DLA in a single jet vs jet energy E for fixed 
Qc = Qo (representing hadrons) and for fixed energy E (Iq = ln(-E'/Qo) = 5) 
but variable resolution Qc (or ycut = [Qc/'^EY) (representing jets) and fac- 
torial moments Eg vs. energy E {¥)■ The moments approach the asymptotic 
limits Fa = 1.33, F4 = 4.62, Fg = 359.7, Fie = 2 x 10^ 



This equation can be solved in terms of Bessel functions |3I 

NiYc) = 2P^Yc + A [h{2P^Yc + X)K^{2pVX) + Ki{2P^Yc + \)h{2l3y/\) 

(18) 

This solution is shown in Fig. HI for constant cut-off Qc = Qo as a full line, 
which represents the well known rise of multiplicity, N ~ exp {2/3 ^/eJA), and 
follows from the asymptotics Iu{z) ~ e^. 

One can also study the variation of jet multiplicity at fixed jet energy E 
with resolution Qc, as discussed previously [16 . This dependence is given 
by the same function (fTH|) if the Durham algorithm is applied; in this case 
the cut-off in the evolution of the parton cascade is defined through the 
transverse momentum > Qc- The corresponding result is represented by 
the dashed line in Fig. ^ One observes a considerably lower multiplicity 
by an order of magnitude, but, for decreasing Qc Qo { Yc ^ ^0) in the 
transition from jet — hadron, the jet curve rises rapidly and reaches the 
curve for hadrons. This behaviour is a consequence of the running coupling, 
which is always smaller for jets, because kx, the argument of a^, is always 
larger. The analytic behaviour in this limit can be derived from (fTHj) using 
the approximation Kq{z) ^ ln(2/2;) for small z with small Qc or large Yc (see 
also t21i) 

N ~ iP^YcY^Hn exp 2/3 v^. (19) 
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This multiphcity diverges logarithmically for Qc ^ A (Ac —>■ 0). This di- 
vergence is a consequence of the diverging coupling in this limit. Since 
<5c > <5o > A, this divergence is outside the physical region, but a strong 
variation remains, as seen in Fig. CI Note that in the case of fixed coupling, 
the moments would only depend on the ratio E/Qc as no other dimensional 
quantities exist; then both curves for jets and hadrons would coincide and 
follow asymptotically a power law ~ (E/Qq)^'^. 

Next we calculate the factorial moments for g > 1. By appropriate dif- 
ferentiation of the evolution equation (fTTjl one finds for g > 1: 




fM = 0. (21) 

where ()2H1 follows from ()12|1 and (jZj), in agreement with (j3)). 

This coupled system of integro-differential equations is solved numerically. 
The integral is computed using the trapezoidal rule with a step size of 6y = 
10~^. As a test, we calculated the multiplicity from the analytic formula and 
compared with the numerical solution of ()16|1 . There was good agreement to 
within 10~^. Reducing the step size by a factor of 10 {6y = 10~^), resulted 
in a change of about 1% in the value of Fiq and a change in the ratio Hr^ by 

0. 1% below in the range considered. 

A selection of results on the normalized moments Eg is also shown in 
Fig. H The moments Eg for g > 1 vanish at threshold = 0. Shortly above 
threshold the moments of higher rank g are suppressed as the multiplicity is 
low here and only a small number of terms contribute to the sum (j3]) ; on the 
other hand, they approach larger asymptotic values. In DLA one finds j2] 

F, = ^y2('^)^^^^^ foTq>l with Fi = l, (22) 
^ m=l ^ / 

1. e. for the first terms 

F_9 F_ 208 2425 _ 2207 

Therefore, with increasing energy more and more moments Eq are needed to 
represent the multiplicity distribution. 

Alternatively, one can study the ratios Hq = Kg/ Eg of cumulant and 
factorial moments. An infesting feature of these ratios is their asymptotic 
decrease with rank g |221 

Hq ~ 1/g^ (24) 
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Figure 2: The function y = Ashio(x) = Arsinh(x/2)/ In(lO) represents data 
with positive and negative signs over a large range of scales, in comparison 
with the common y = log^^Q x. 

This behaviour can be derived easily in leading order of 70 ~ ^/a^ under 
the assumption of constant Kg and Fg from the equation kg = 7o/(j, which 
follows from ((7j) and (fTTj) using (fTTjl . 

On the other hand, these ratios are not very convenient near threshold. 
First, we can derive for the cumulant moments from and 1)211) 

Kg = {-iy-\q-iy. forre = 0. (25) 

Then, because = for g > 1, the ratios Hg diverge in the limit Yc ^ 
with alternating signs. Note that the observables Fg,Kg, Hg each provide 
a full description of the multiplicity distribution. At high energies Kg and 
Hg are convenient, whereas at low energies Fg is more convenient since the 
respective higher order terms are suppressed. 

Before we come to a discussion of the results, we introduce a convenient 
way to present data with a large range of scales and alternating signs. It is 
an extension of the usual log plot for positive numbers in which a positive 
quantity y is represented in log scale y = 10^, x = log^Q?/. For a quantity 
with either sign we write 

y = 10^ - 10"^ = 2Sinh(x In 10) (26) 
X := Ashio(y) = (1/ In 10)Arsinh(|//2) (27) 

where Arsinh(2) = ln(2; + vl+^). For large positive or negative numbers 
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Figure 3: Moment ratios Hq as function of rank q in DLA for different 
energies Y for fixed Qc = Qo (representing hadrons). The quantity plotted is 
Ashio(ifgXlO^) with function With increasing Y the moments approach 
the asymptotic behaviour ifg ~ 



one obtains Ashio(±?/) ±log]^o(?/) For convenience we show in Fig. |21the 
comparison of both functions. 

First, we study the evolution of the ratios, Hq, with energy Y = \n{E/Qo) 
and this is shown in Fig. El for a few values of Y. One can see the oscillations 
with large amplitude near threshold (small Y) which follow from ()25p. With 
increasing energy the oscillations continue with the same oscillation length 
but with decreasing amplitude. The moments of even rank q finally change 
sign and the oscillations disappear at Y = 1 where they are already close to 
the asymptotic limit (jUJ. Also one observes the rise of Hq with q at small 

Y and the decrease of Hq with q at large Y. 

Next we present in Fig. H] the evolution of the ratios Hq with energy 

Y which diverge at threshold and approach the asymptotic result (j^^ from 
below. At the same time we show the variation of the multiplicity moments 
for jets at variable resolution at fixed energy Y = 5 (corresponding to LEP 
energies at Qq ~ 0.3 GeV). Similarly to the case of multiplicity, there is a 
large difference between both dependences, which reflects again the role of 
the running coupling. The rise of the jet moments is delayed but they reach 
the hadron moments as Qc — »■ Qo at the nominal energy oiY = 5. 
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Figure 4: Ratios of moments Hq in DLA vs energy for fixed Qc = Qo 
(representing hadrons) and for fixed energy Y = ln{E/QQ) = 5 but variable 
Qc representing jets at variable resolution ycut = (Qc/'^E)'^. The asymptotic 
limits Hq ^ are shown as well. 



Finally we remark that the comparison with e^e~ results would require 
quark jets in two hemispheres. This can be obtained |S| by replacing the 
generating function as 

Z ^ Z^^^/^-* = (28) 
We do not go into this minor difference in the present qualitative discussion. 



3.2 Modified Leading Logarithmic Approximation 

In this approximation, the full DGLAP splitting functions |37l38j are in- 
cluded, with energy conservation taken into account. The restriction to 1- 
loop results remains. MLLA takes into account the next-to-leading order 
terms in the ^fa's expansion of the exponent, InA^^. In the present analysis, 
we neglect the qq pair production (even though it contributes in MLLA) for 
simplicity - it could be taken into account without difficulty (see, for exam- 
ple |l(ipi7j ). The evolution equation for the generating function Z{YciU) in 
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gluodynamics reads [T3j 



dYr, 



J Zn 



x{Z(rc + lnz,u)Z(Fc + ln(l - z),u) - Z(y;,w)}(29) 
Z(0,m) = u (30) 



where 



~ 1t{ ~ 
Zc = e'^" , asikx) = = , kx = mm{z, 1 — z)E (31) 



b \n{kT/A) 

and the sphtting function Pgg{z) is given by^ 



, , 1 - z z 
z(l -z) + + 



1-z 



(32) 



For the symmetric kernel one can replace in the integral ^Pgg{z) by Pgg'^{z) = 
(1 — z)Pgg which has no pole at ^ = 1 and behaves for z ^ like P^^^ ~ 6/z. 
Differentiation by u leads again to evolution equations for the multiplicity 



N'g{Y^)= dz^^P-y{z)x 

X {NgiY, + Inz) + NgiY, + ln(l - z)) - N,{Y,)} (33) 



and the unnormalized factorial moments 



/;ra= l''^^dz^^p-y{z)x 



X {/,(F, + In^) + /,(y, + ln(l - z)) - /,(F,) + 

+ E (^)^"(^- + l^^)V™)(^c + ln(l-;2))} (34) 

m=l ^ 

At threshold for multiparticle production, we find the following conditions 

q = l: Ng = 1 ioi E < 2Q^, q > 1 : fg = OioiE< qQ^. (35) 

The condition for g = 1 follows from the initial condition Ng = 1 and N'g = 
in ()33|) below the threshold for particle emission at E' < 2Qc. 

^The P function "371 is related to the $ function ^ by <^ab = 2Pafc 
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The constraint for the higher moments, fg for E < qQc, corresponds to the 
kinematic condition for producing q particles. We can derive the condition 
for g > 1 by induction. First, consider the case of g = 2 in ()34|1 . We know 
f^(Yc) = if Fc < ln2 (equivalently, ii E < 2Qc) from fi=N = lm this 
region. Now let us consider the situation for q > 2, taking for g — 1 
as given. Consider first the evolution of fq{Yc) in ()34|) with Yf. starting from 
fq = a.t Yc < \n2 where N = 1. As long as the sum in ()34|) vanishes we 
just have the evolution fq{Yc) = fq{Yc) x I(Yc) with a known integral / and 
therefore fq = 0. 

Next we show that for energies E < qQc (or Yc < Inq) this sum indeed 
vanishes. First, note that the argument hi[zE/Qc) of the moment fm varies 
from at the lower limit to ln(g — 1) at the upper limit while z varies from 
to 1 — Zc and the argument of f[q-m) decreases from ln(g — 1) to 0. In the 
sum, exactly one of the factors - fm at energy E' = zE or fq-m at energy 
E" = {1 — z)E - vanishes because of condition ()35j) for m < q ot q — m < q 

E'<mQ, E">{q-m)Q, U = > (36) 

E'>mQ, E"<{q-m)Q, = > (37) 

Now, the perceptive reader may notice that for m = 1, fm = Ng cannot be 
0. This is no problem, however, since fq-i = in the whole range. Thus, we 
have proven the validity of for any q. 

The MLLA evolution equation reduces to the DLA equation (fTTj) 
if the integrand is taken in small z approximation, i.e. the ln(l — z) term 
is neglected and Pgg{z) ~ This limit should be achieved for very high 
energies, because ^^(/ct) remains energy independent and gives weight to ever 
smaller values of z under the integral. Indeed, the high energy results for the 
moments in terms of an ^/a^ expansion |ll|22j yields the DLA solutions in 
the asymptotic limit. 

Full solutions of the evolution equations (j33|) . (j3l|l satisfying the threshold 
condition, analogous to (jlHj) in DLA are not known. If the integrand is 
simplified assuming a ^/aj expansion one obtains expressions with Bessel 
functions similar to (fTHj) but with non-integer index ^9 39j which remain finite 
for Qo A. 

Here we solve the system of equations again numerically using 

the trapezoidal approximation with step size 10~^ for the integral (for control 
also with 10~^), alternatively the 3-point Simpson rule has been applied with 
similar results. Note that the evolution variable Y^ ~ InE now not only 
appears in the integral bound but also in the integrand which was not the 
case in the DLA. 

In Fig. El we show the numerical results for the multiplicity in a single 
jet as function of Yc for hadrons (fixed Qc = Qo according to LPHD) and for 
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Y = ln(E/QO) 

Figure 5: Multiplicity in DLA in single jet vs energy for fixed Qc = 
Qo (representing hadrons) and for fixed energy Yq = ln(-E/Qo) = 5 but 
variable Qc representing jets at variable resolution Qc (or i/cut = {Qc/'^EY) 
and factorial moments Fq vs. energy Yq. 

jets at fixed energy E (LEP energy) but variable resolution. The results look 
similar to Fig. ^in DLA, but the absolute size is now much reduced because 
of the energy conservation constraint. 

Also shown are some factorial moments, Fq. Again, the moments are 
smaller in size and, therefore, the approach to the asymptotic values (j23|) is 
further delayed. A distinctive difference from DLA are the shifted threshold 
energies for the higher moments according to 

The factorial moments Fq show a striking result at intermediate energies: 
they all cross at one point {E ^ 20 GeV) with = 1, i.e. with a Poissonian 
distribution 

Poissonian transition point: Fq = 1 at Yp 4.7. (38) 

This feature does not appear in DLA, and is apparently related to the delayed 
thresholds. For the considered moments with g < 16 the crossing at = 1 
occurs within an interval of 5Y ^ 0.05 for both step sizes 10~^ and 10~^ 
of our numerical computation, so the distribution is close to a Poissonian. 
We note though, that the higher moments Fq with g > 32 have a threshold 
above the transition point Yp so the Poissonian cannot be an exact solution. 
Nevertheless, it is remarkable that besides the threshold and asymptotic 
regime, there is an intermediate energy with a very simple behaviour. The 
Poissonian transition point follows apparently from the evolution equation - 
a novel property for which we have no analytical explanation yet. 

Next we consider the moment ratios Hq for hadrons {Qc = Qo) from low 
to high energies in Fig. IHl At the low energies {Y < 4) we observe again 
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Figure 6: Moment ratios Hq as function of rank q in MLLA for different 
energies Y for fixed Qc = Qo (representing hadrons). With increasing Y the 
moments approach the asymptotic behaviour Hg ~ l/q'^- 

the pattern of rapid oscillations from one order q to the next as in DLA 
reflecting the threshold behaviour. Beyond the Poissonian transition point 
[Y > 6) we enter a new regime with a larger oscillation length. The first 
minimum is slowly rising from g = 5aty = 6tog = 7aty = 10. At 
higher q a maximum follows. This is expected from the analytical results 
|13|22j . Whereas the oscillation amplitude of Hq rises with q at low energies, 
it decreases at the high energies. However, it is always much below the 
asymptotic DLA limit for large q, only for small q below the first minimum 
one can see a convergence to this limit in the considered energy range. 

The evolution of the Hq moments with Y^ is shown in Fig. [7|for hadrons 
and jets. In comparison to the DLA results of Fig. ^ we see at first a 
qualitatively similar behaviour, but again the approach to the asymptotic 
limit is not visible for the higher moments in the considered energy range. 

More details are shown for hadrons in Fig. |H1 New features emerge 
at high energies (above Y ^ Yp ^ 4). For the higher moments {q > 4) 
secondary minima and/or maxima appear with increasing energy at the level 
of O{10~'^). These secondary oscillations lead to a delay of the onset of the 
asymptotic behaviour. Below Yc = 20 {E < 10^ GeV) the moments with 
q >6 are below 1/10 of their DLA asymptotic value. 
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Figure 7: Ratio of moments Hq in MLLA vs energy for fixed Qc = Qo 
(representing hadrons) and for fixed energy Y = ln{E/QQ) = 5 but variable 
Qc representing jets at variable resolution ycut = (Qc/'^E)'^. The asymptotic 
limits Hq ^ are shown as well. 
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resenting hadrons) as in Fig. [3 The asymptotic limits of Hq are above the 
upper bound of the plot. 
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The Poissonian transition point corresponds to Hq = for all g > 1. In- 
deed the moments show this behaviour in good approximation at the above 
Yp. Beyond this point the positive short range correlations lead to a positive 
H2, whereas below it the energy conservation constraints lead to a negative 
cumulant. In general, the even Hq ratios change sign because of the neg- 
ative value at threshold and the positive asymptotic limit. In hadron 
phenomenology, the positive correlations H2 > reflect the resonance pro- 
duction jlOj, in our calculation it is the gluon bremsstrahlung with small 
angle (above Qo cut-off) which lead to the short range correlations, be it 
between "hadrons" (partons at scale Qo) or between (mini) jets. 

As common features of DLA and MLLA, we note the splitting between 
hadron and jet moments at the same which is a direct consequence of the 
running coupling. Hadron and jet multiplicities coincide for Qc ^ Qo- 

4 Monte Carlo simulation of parton cascade 

After the numerical solutions of the DLA and MLLA evolution equations 
for single jets, we finally apply an MC event generator (ARIADNE [32\ ) 
at the parton level based on the same procedures as the above evolution 
equations: perturbative QCD evolution with the coupling asikx) terminated 
by a transverse momentum, kx, cut-off and arbitrary parameters A and Qq > 
A. The MC involves the coupled evolution of quarks and gluons in the 
cascade, the inclusion of large angle radiation, the full first-order matrix 
element for e~^e~ — > qqg, and exact energy- momentum conservation. These 
features lead to a higher accuracy than our MLLA calculation.^ On the 
other hand, a simplification of the Monte Carlo is the exclusive use of 1- 
loop calculations. This simplification could result in deviations when very 
different energy scales are compared. 

In our adoption of the MC program, we set all quark masses to zero 
but kept the masses for heavy quarks in the calculation of the number of 
flavours nj for «<, at a given dipole mass. The results depend only on the 
two adjustable parameters A and Qq; alternatively, we use the parameter 
A = ln((5o/A). These parameters have been determined from a fit to multi- 
plicities. We refer to this modified Monte Carlo as "ARIADNE-D," where 
"D" stands for "Duality" . 

The jet multiplicities are obtained from the final state partons using the 

^In MLLA, the inclusion of quark jets is possible - also, the first-order matrix element 
has been included in the calculation of multiplicities 16 , however, energy- momentum 
conservation is fully realized in the MC approach and large angle emission is more easily 
accessible. 
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Durham algorithm, just as in the experimental analysis with (jHJ and ©• 
Hadron multiplicities are related to the parton multiplicities, generated at 
scale Qo, according to the LPHD prescription. 

This model will be compared with the experimental data. First we con- 
sider the mean hadron and jet multiplicities A^, as shown in Fig. IHl The jet 
multiplicities for energies Q = 35, 91 and 189 GeV are obtained as functions 
of the cut-off Qc or of i/cut = {Qc/QY ^ in the figure they are plotted as a 
function of the variable 

with the additional scale Qq. The MC parton multiplicity Npart is shown 
at the cut-off Qc = Qo and should be compared at this scale with the to- 
tal hadron multiplicity Ntot according to the LPHD prescription, in general 
Ntot = K X Npart (where i^' is a constant); the observed charged multiplic- 
ity is given by N^h = fchNtot with the charged particle fraction fch- These 
factors also depend on whether weak and electromagnetic decays of hadrons 
are included {K^ ., A, . . .). The factor Kfch will be determined from the fit. A 
simple choice is = 1 as for jets and fch = "2/3 a.s obtained previously |16].^ 
For jets we compare parton and hadron results at the same cut-off and with 
K = 1. Note that the choice of variable does not matter for the comparison 
of jet properties (both MC and experimental data are obtained by the same 
procedure); but, using the variable Yes in (jHTIjl . has the advantage that the 
jet data approach the hadron data for ycut — at a finite value of Yes in the 
figure according to In the region Qc = I — 2 GeV the dependence on 
Qo disappears and we are in the perturbative QCD regime, governed by the 
single parameter A |52j . 

In order to determine the parameters of the model we compared the 
MC results with the multiplicity data at LEP-1. We noted that we could 
not get good agreement with the jet data from OPAL j3J over the entire 
kinematic range of 10~^ < Ucut < 0.5. Since we were interested in studying 
the transition between jets and hadrons, we chose our parameters so that they 
gave a particularly good description of the low ycut regime and the hadron 
multiplicity as well as the very large ycut region. In Table Q we illustrate 
the dependence of the hadron and jet multiplicity on the parameters. We 
require a good fit of the jet multiplicity of A^ = 14.1 ± 0.4 for ycut = 2 x 10~^ 
jjl] corresponding to Qc = 0.4 GeV. At this scale the particles from non- 
hadronic decays are recombined into the primary hadrons and this ambiguity 

''if proportionality of energy spectra of partons and hadrons according to LPHD is 
required in the fuU energy range, then the normalization has to be the same (K — 1) 
because of energy conservation. 
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Figure 9: Multiplicity N of hadrons (taken as Net x 1.25 using data 
|41 l4'JI4,'AI44l45l4fjl4VI48l4flHl] and LEP averages from [22]) and multiplicity 
of jets |51|28j in e~^e~ annihilation together with our Monte Carlo calculations 
as function of Yes, see flS^ (F^s ~ ln(Q/Qc))- 



is removed. Also we considered the mean charged multiplicity at LEP-1 of 
Nch = 21.2 ±0.3 j22j (this number includes particles from and A decays). 
Solutions can be found close to the simple case fchK = 2/3. We do not 
require this condition though and choose the parameters 

A = 400 MeV, A = 0.01, (KU)"^ = 1.25; (40) 

they are not very different from the previous study for hadrons alone (A = 200 
MeV, A = 0.015 and {Kfch)~^ = 1-5 |53j) which optimized the description of 
the hadronic correlations alone. 

With our choice ()40p a reasonable description of jet multiplicities for small 
and large ycut was obtained, not only at LEP-1, but also at the higher and 
lower energies of 35 and 189 GeV, as can be seen in Fig. IHl Furthermore 
we get a good description of the trend of hadron multiplicities as function of 
cms energy in the full range Q = 3 . . . 200 GeV. The Monte Carlo was run 
with 250,000 events at each energy. 

As in the analytic calculations, the multiplicity in the MC approach ap- 
plied here would diverge for Qo — > A. Now, Fig. IHl shows that with the same 
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Table 1: Monte Carlo results on hadron and jet multiplicities for given values 
of A (rows) and A (columns). The first half of the table consists of data with 
Vcut = (hadrons), and the second half consists of data with i/cut = 2 x 10~^. 



A(MeV): 


A = 0.001 


A = 0.01 


A = 0.015 


A = 0.05 


250 


35.8 


31.7 


30.6 


26.5 


300 


33.3 


29.5 


28.4 


24.8 


400 


29.6 


26.3 


25.4 


22.0 


500 


27.0 


24.0 


23.2 


20.1 


250 


12.4 


12.2 


12.1 


11.8 


300 


13.0 


12.8 


12.7 


12.4 


400 


13.9 


13.6 


13.5 


13.0 


500 


14.4 


14.0 


13.9 


13.2 



cut-off Qo taken at all energies the correct energy dependence of the hadron 
multiplicity is obtained. 

At every cms energy in Fig. IHl the MC fails to duplicate the intermediate 
Ucut region with MC results falling above the experimental data. We find 
that the structure we are unable to reproduce occurs at the fixed b mass 
scale Qc ~ 5 GeV at all energies. In our Monte Carlo we use only light 
quarks so as to avoid complicated decay processes. It is plausible to assume, 
then, that the 'hump' structure at 5 GeV comes from b quark production. 
This conjecture is further supported by the enhancement of the effect seen at 
the Z energy of 91 GeV as expected from the enhanced b quark production 
by neutral current interactions of the Z. 

The results in Fig. IHl look similar to the analytic solutions in Figs. C] 
and El and demonstrate again the importance of the running coupling which 
leads to the scale (Q) dependence of multiplicity at fixed Fes or Q/Qc] in the 
threshold region (low Y^s) the difference between hadron and jet multiplicities 
amounts to a an order of magnitude. 

Next, we turn to the energy evolution of the Hg moments for hadrons 
(i.e. partons at scale Qo) from low to high energies (10-1000 GeV). This 
evolution is displayed in Fig. Eland looks similar to the result from MLLA 
in Fig. ini at the low energy, below the Poissonian point, there are the 
rapid oscillations reflecting the threshold behaviour, above that point we 
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Figure 10: Ratios Hg for a sequence of cms energies Q in e~^e annihilation 
as obtained from the Monte Carlo calculation, plotted as in Fig. El 



observe the oscillations with increasing length for increasing energy. The 
first minimum occurs at qmin ~ 5 at 90 GeV and increases to g^m ~ 7 at 
1000 GeV. This increase is similar but a bit stronger than in the asymptotic 
formula (jHl), which would predict an increase by one unit of qmin instead of 
two. These calculations have been carried out with a sample of 4M events. 
The statistical errors are determined from the fluctuations of 4 sub-samples 
and are found to be significantly smaller than the oscillation amplitudes at 
their maximum. The numerical results are also given in Table 121 

In Fig. ^2 we compare the ratios, Hg, obtained from our calculations 
with the measurements by L3, both for hadrons and for jets with variable 
(Durham) cut-off Qc- For hadrons we assume that the multiplicity distribu- 
tion of charged particles has the same normalized moments as the distribution 
of all particles. This is the case if the same conversion factor fch applies for 
all multiplicities n. 

For hadrons, a measurement of the full multiplicity distribution j22l| which 
extends up to rich = 62 and one with truncation [2H] at the multiplicity 
rich = 48 are presented. The truncation can be seen to considerably increase 
the oscillation amplitude. We compare with the Monte Carlo results again 
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Table 2: Ratios of multiplicity moments, Hq, for different cms energies, Q, 
in e^e~ anihilation obtained from ARIADNE-D parton MC (the statistical 
errors of the MC results are generally smaller than the last digit given) 



Q 


10 GeV 


35 GeV 


91 GeV 


190 GeV 


1000 GeV 


2 


-0.0565 


0.00808 


0.039 


0.0569 


0.0864 


3 


0.0129 


0.00495 


0.0074 


0.0104 


0.0176 


4 


-0.0047 


-0.00064 


0.00040 


0.00126 


0.00356 


5 


0.0022 


-0.00015 


-0.00022 


-0.00014 


0.00040 


6 


-0.0013 


-0.000051 


-0.00012 


-0.00020 


-0.00019 


7 


0.0010 


0.000028 


-0.000017 


-0.00008 


-0.00021 


8 


-0.0010 


0.000015 


0.000013 


-0.000006 


-0.00014 


9 


0.0011 


-0.0000005 


0.000006 


0.000015 


-0.00008 


10 


-0.0015 


-0.0000069 


-0.000003 


0.000009 


-0.00004 


11 


0.0022 


-0.0000041 


-0.000005 


-0.0000009 


-0.000005 


12 


-0.0038 


0.0000032 


-0.000001 


-0.000005 


0.000027 


13 


0.0071 


0.0000036 


0.000003 


-0.000003 


0.000047 


14 


-0.013 


0.0000014 


0.000004 


0.000001 


0.000045 


15 


0.026 


-0.0000048 


0.000001 


0.000005 


0.000019 


16 


-0.035 


-0.0000029 


-0.000003 


0.000003 


-0.000018 


17 


-0.09 


0.0000008 


-0.000005 


-0.000001 


-0.00005 


18 


1.7 


0.0000093 


-0.000002 


-0.000005 


-0.00005 


19 


-19 


0.0000018 


0.000003 


-0.000005 


-0.00003 


20 


250 


-0.0000093 


0.000007 


0.000001 


0.00002 



with and without truncation. The truncation was at multiplicity nt = 63 to 
obtain the same event fraction of 0.005% for multiplicities n > as in the 
experimental data. Our results reproduce the data rather well but for the 
truncated moments the high q oscillations are a bit weaker. 

We mention here another subtlety of the experimental analysis. The 
L3 data in Fig. [TT] refer to data without and A decays, whereas the 



21 




-0.8 



xlO 
0.05 



-0.05 
-0.1 

-2 

xlO 
0.05 


0.05 
0.1 
0.15 





4 6 8 


10 12 

rank q 




Qc=100 


MeV 




\ /* 0.04 
* 0.01 


• • 








-0.02 


1 1 1 




2 3 4 



14 16 1i 



0.025 



Qc=400 MeV 




-0.025 



-0.05 - 



0.5 - 



- 




Qc=3GeV 



2^ 3^ 4 I 

5 10 15 5 10 

rank q rank q 

Figure 11: (a) Ratio Hg of multiplicity moments obtained by L3 Collab- 
oration |28|29j with truncation of multiplicity at rich = 48 compared with 
ARIADNE-D MC (truncation at = 65) and both without truncation; (b) 
Ratios Hg for jets at different cut-off Qc [2H] compared with our MC results. 
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Figure 12: Ratio of moments Hq for hadrons IMJ and jets |28|29j in com- 
parison with ARIADNE-D Monte Carlo. 



multiphcities in Fig. El include them. However, the Hq moments do not 
depend on this difference, except to within a few per cent. The same is true 
for the truncated moments, as long as the removed event fraction is the same 
|29j . This appears plausible from the point of view of KNO scaling. 

We note that the moments presented by the SLD Collaboration |2Z] have 
been determined from the multiplicity range 6 < rich < 54 and they look 
very similar to the truncated L3 moments. We stress that from our point of 
view the untruncated moments are of primary relevance; only these obey the 
evolution equations discussed above and have simple asymptotic properties. 



23 




Beyond the first minimum they are very small at LEP energies, of the order 
of 10"^ 

Finally, in the same Fig. ^2 we study the evolution of the Hq moments 
with cut-off Qc or Ucut = {Qc/Q)"^- The MC calculations reproduce well the 
change from the rapid oscillations with large amplitude at large cut-off Qc/ Q 
to the oscillations with larger oscillation length and smaller amplitudes in the 
small Qc/Q region close to the fully resolved hadronic final state. 

In the subsequent Figs. IT^HT^ we show these ratios, Hg, separately as 
function of Fes, i-e. of Q/Qc- Included are further data at intermediate i/cut 
values from [23 together with the corresponding MC predictions and also the 
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Figure 14: Ratio of moments Hq for hadrons and jets as in Fig. El 

predictions for the hadronic final state, i.e. at Qc = Qo- In these calculations 
2.5 M events have been generated in the MC. For all orders g, the calculation 
follows the experimental data closely. Again, there is the same qualitative 
difference between the hadron and jet data as observed before for the DLA 
and MLLA calculations in Figs. |3]and[7| 

For the ratio, we also compare with results we deduced from TASSO 
moments 54 . The data clearly confirm the splitting of the moments at 
different scales, as expected from the running coupling. At small Q/Qc-, both 
for hadrons and for jets, the distribution is narrower than Poisson {H2 < 0), 
while at higher values of this ratio the distribution gets broader {H2 > 0). 
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The positive Hg is expected from the short range correlations, either from 
gluon bremsstrahlung at small angles for partons or from resonance decays 
for hadrons. At the point with H2 = the higher moments are also close 
to zero corresponding to a Poissonian transition point, both for jets and 
hadrons: 



It is interesting to note that the Poissonian distribution has actually been 
successfully fit to their data at 29 GeV by the HRS Collaboration [45J. 

We observe again a close similarity between MC and MLLA results, al- 
though the deviations from a Poissonian {H^ > at minimum) appear a bit 
larger in the MC. 

5 Conclusions 

We have studied the multiplicity moments, in particular the Hg ratios for 
the QCD parton cascade in three different approximations: DLA, MLLA 
and parton Monte Carlo. The DLA predicts the very asymptotic behaviour. 
The asymptotic MLLA corrections predict the existence of Hg oscillations, 
but do not explain any i/cut dependence. The more precise calculations of the 
full MLLA solution of the evolution equations and the Monte Carlo method, 
which include the correct threshold behaviour, show at the presently acces- 
sible energies the following new features of the multiplicity moments: 

a) There is a particular energy scale, Yp, for hadrons or a particular ycut 
for jets where the multiplicity distribution is approximately Poissonian 
{Hg ^0, Fg^lfoT q> 1). 

b) In a region of smaller energies there is a rapid oscillation of Hg, reflect- 
ing the threshold behaviour; at higher energies the oscillation length 
increases with energy as predicted qualitatively by asymptotic MLLA, 
but the oscillation amplitude is much smaller. In fact, for hadrons at 
LEP energies we predict \Hg\ < 10~^ beyond the first minimum at 
q ^ 5. Also, there are secondary extrema of Hg as a function of energy 
for a given rank q. The approach to the asymptotic DLA region is very 
slow if visible at all. 

c) For fixed resolution parameter ycut the moments of jets and hadrons 
are rather different as a consequence of the running coupling. 



Poissonian transition point: 




Q = 91, Qc ^ 0.3 GeV jets. 



hadrons 



(41) 
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We conclude that the perturbative approach yields good results not only 
for single particle phenomena, such as single inclusive spectra and mean mul- 
tiplicities, as envisaged originally in the LPHD approach, but also for corre- 
lations of high order. Above the Poissonian transition point, the short range 
correlations initiated by resonance decays for hadrons and gluon Brems- 
strahlung for partons become increasingly important and lead to the broad- 
ening of the multiplicity distribution of both hadrons and jets {H2 > 0). So 
there is a duality also between hadronic and partonic correlations. 

Moreover, the absolute normalization, originally considered as free param- 
eter, can be studied as function of resolution i/cut and a unified description of 
hadronic and jet phenomena is possible, using the variable Yes in (jH^ with 
the shifted t/cut prescription, both for multiplicities and correlations. 

Both the predicted energy evolution of hadronic correlations and the ycut 
dependence of the jet correlations follow the experimental results closely, with 
the exception of jets near the b quark threshold. In the jet regime Qc > 1 — 2 
GeV, the effect of the hadronic scale Qq becomes negligible and we are in the 
domain of perturbative QCD with the single parameter A only; Q^. is then 
an external parameter to be chosen by the experimenter. At smaller scales, 
the hadronic scale Qo becomes important, and the results become model 
dependent. The transition region from jets to hadrons for i/cut —>■ 0, which 
shows strong variations of all moments, is very well reproduced by the parton 
MC with low cut-off > Qq > A. In this region the coupling becomes large, 
i.e. cts > 1. An extension of the calculation into this kinematic region is not 
justified a priori, but, as seen from the analytical results, the convergence of 
the perturbation theory is not in danger, and the solutions can be obtained 
from the evolution equations or from the MC. This successful description can 
be viewed at least as a good effective parametrization of the soft transition 
region. 

With increasing order, accurate calculations are required. So far, this is 
only accessible by numerical methods. Further improvements are possible by 
including heavy quarks in the parton evolution and by including the 2-loop 
results. 

The good description of the data by the MC model with normalization 
K ^ 1 implies that in the considered applications, one parton counts about 
one hadron and the hadron final state is well represented by a parton final 
state of the same multiplicity at resolution scale Qo > A. 
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